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Abstract 

O 
O 

The quark-antiquark potential is constructed by making use of a new analytic 
running coupling in QCD. This running coupling arises under "analytization" 
of the renormalization group equation. The rising behavior of the quark- 
■ antiquark potential at large distances, which provides the quark confinement, 

is shown explicitly. At small distances, the standard behavior of this potential 
originating in the QCD asymptotic freedom is revealed. The higher loop 
CN ' corrections and the scheme dependence of the approach are briefly discussed. 
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I. INTRODUCTION 



The description of quark dynamics inside hadrons remains an actual problem of elemen- 
tary particle theory The asymptotic freedom in quantum chromodynamics (QCD) enables 
one to investigate the quark interactions at small distances by making use of standard pertur- 
bation theory. The quark dynamics at large distances (the confinement region) lies beyond 
such calculations. For this purpose other approaches are used: phenomenological potential 
models string models 0, bags models [[|, lattice calculations the explicit account of 
nontrivial QCD vacuum structure ||, and variational perturbation theory ||. 

Recently Shirkov and Solovtsov proposed a new analytic approach to QCD |7| . Its basic 
idea is the explicit imposition of the causality condition, which implies the requirement of the 
analyticity in the Q 2 variable for the relevant physical quantities. The essential merits of this 
approach are the following: absence of unphysical singularities at any loop level, stability 
in the infrared (IR) region, stability with respect to loop corrections, and extremely weak 
scheme dependence. The analytic approach has been applied successfully to such problems 
as the r lepton decays, e + e~-annihilation into hadrons, sum rules (see || and references 
therein) . 

In Refs. P,[ToH the analytic approach has been employed to the solution of the renormal- 
ization group (RG) equation. The analyticity requirement was imposed on the RG equation 
itself, before deriving its solution. Solving the RG equation, "analytized" (i.e., requiring 
analyticity) in the above-mentioned way, one gets, at one-loop level, a new analytic running 
coupling, which possesses practically the same appealing features as the Shirkov-Solovtsov 
running coupling does. An essential distinction, that will play a crucial role in the present 
paper, is the IR singularity of the new analytic running coupling at the point q 2 = 0. It 
should be stressed here that such a behavior of the invariant charge is in a complete agree- 
ment with the Schwinger- Dyson equations, and, as it will be demonstrated in Sec. Ill, 
provides the quark confinement (see Sec. II for the details). 

In this paper we shall adhere to the model P,|TT| of obtaining the quark- ant iquark (qq) 
potential by the Fourier transformation of the running coupling. However, the perturbative 
running coupling a s (q 2 ) does not enable one to obtain the rising qq potential without invoking 
additional assumptions [|TTj| . 

The objective of this paper is to construct the quark- ant iquark potential by making use of 
the new analytic running coupling. This potential proves to be rising at large distances (i.e., 
providing the quark confinement) and, at the same time, it incorporates the asymptotic 
freedom at small distances. It is essential that for obtaining this potential no additional 
assumptions, lying beyond the standard RG method in the quantum field theory and the 
analyticity requirement, will be used. 

The layout of the paper is as follows. In Sec. II the derivation of the new analytic running 
coupling is presented and its properties are briefly discussed. In Sec. Ill the quark- ant iquark 
potential, generated by the new analytic running coupling, is derived by making use of the 
Fourier transformation. Further, the asymptotic behavior of the potential at large and small 
distances is investigated. In Sec. IV the higher loop corrections and the scheme dependence 
of the potential are discussed briefly. For practical purposes, a simple approximate formula 
for the potential is proposed which interpolates its infrared and ultraviolet asymptotics. 
This formula is compared with the phenomenological Cornell potential. Proceeding from 
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this, an estimation of the QCD parameter A is obtained. In the Conclusion (Sec. V) the 
obtained results are formulated in a compact way, and the further studies in this approach 
are outlined. 



II. A NEW MODEL FOR THE QCD ANALYTIC RUNNING COUPLING 

In the analytic approach to QCD, proposed by Shirkov and Solovtsov 0, the basic idea 
is the explicit imposing of the causality condition, which implies the requirement of the ana- 
lyticity in the Q 2 variable for the relevant physical quantities. Later this idea was applied to 
the "analytization" (i.e., the procedure of analyticity requirement) of the perturbative series 
when calculating the QCD observables 0. The results turned out to be quite encouraging. 
As was mentioned in the Introduction, the analytization of the perturbative series leads to 
the elimination of the unphysical singularities, to the higher loop correction stability and 
to a weak scheme dependence. However, the Q 2 -evolution of some QCD observables (for 
instance, the structure function moments) is intimately tied with the solution of the RG 
equation. Our task here is to involve the analytization procedure into the RG formalism 
more profoundly. 

Let us consider the RG equation of a quite general form for a quantity A(Q 2 ) (it may 
be, for example, the gluon propagator, or the structure function moment). At the one- loop 
level this equation reads 

dlnA(Q 2 ) 



rflnQ 2 



7&f ) (Q 2 ), (1) 



where 7 is the corresponding anomalous dimension (the negative noninteger number in the 
general case), a^(Q 2 ) = l/\n(Q 2 /A 2 ) is the one-loop perturbative running coupling. The 
solution of Eq. ([[]) can be written in the form 



A(Q 2 ) = A(Q, 



a 



(1) 



(Ql) 



a { s\Q 2 



(2) 



From here it follows immediately, that this solution has unphysical singularities in the physi- 
cal region Q 2 > 0. However, in many interesting cases mentioned above, the quantity A(Q 2 ) 
must have correct analytic properties in the Q 2 variable (namely, there is the only cutoff 
Q 2 < 0). One can demonstrate this proceeding from the first principles. So, for the gluon 
propagator this assertion follows from the causal Kallen-Lehmann representation (see, e.g., 
fl2fl)) and for the structure function moments this is a consequence of the Deser-Gilbert- 
Sudarshan integral representation^] (see, e.g., [T3]]). Thus, we come to a contradiction. 



The point, which is crucial to our consideration, is the following. The RG equation in 
the form ([!]) involves, in fact, a contradiction. The left-hand side of this equation has no 
unphysical singularities in the Q 2 > region, while its right-hand side has pole-type singu- 
larity at the point Q 2 = A 2 . The account of the higher loop contributions just introduces 



In the most general case this follows from the Jost-Lehmann-Dyson integral representation for 
the structure function, but the detailed discussion of this point is beyond the scope of this paper. 
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the additional unphysical singularities of the cut type in the physical region Q 2 > and 
hence does not solve the problem. 

In order to avoid this contradiction, we propose to use the following method P,|10 



Before solving the RG equation ([![) one should analytize its right-hand side as a whole. This 
prescription leads to the analytized RG equation, which, at the one-loop level, takes the 
form 

where a { ^(Q 2 ) = 1 / ln(Q 2 / A 2 ) + 1/(1 - Q 2 /A 2 ) is the perturbative running coupling analy- 
tized by making use of the Shirkov-Solovtsov prescription 0.0 The solution of Eq. @ can 
be presented in the form 



A(Q 2 ) = A(Q 2 ) 



N 5ii } (g 2 )j 



(4) 



where 



N ain(<5 2 ) = -J—1 * = (5) 
z m z \ z 



Comparing the solution (|J) with Eq. (|J) one infers that N a^(Q 2 ) should be treated as a 
new one-loop analytic running coupling. Really, it possesses the same properties, as the 
one-loop running coupling analytized through the Shirkov-Solovtsov procedure. Namely, 
the new running coupling has the standard asymptotic behavior at z — > oo and it has no 
unphysical singularities in the Q 2 > region. The latter follows directly from the causal 
representation of the Kallen-Lehmann type, that holds for N a^(Q 2 ): 

N s£W)=r^^, N p(-)=(i+ 1 ) r ^ i -^. (6) 

Jo u + z V cry In a + tc 2 

The distinctive feature of the new running coupling, which will play the crucial role in 
the framework of our consideration, is its singularity at the point Q 2 = 0. It is worth noting 
that such a behavior of the invariant charge is in a complete agreement with the Schwinger- 
Dyson equations (see discussion in Ref. 0]), and, as it will be demonstrated in the next 
section, provides the quark confinement. 

Summarizing all stated above, we propose the following model for the analytic running 
coupling. We define the new analytic running coupling N a an (Q 2 ) as the solution of the 
analytized RG equation at the respective loop level. Here one has to choose the anomalous 
dimensions in such a way that the solution of the standard RG equation is the perturbative 
running coupling at the loop level considered.^ Thus, at the one-loop level, the new analytic 
running coupling has the form PJTD[| 



2 It is worth noting here that there is no consistent way for analytizing the RG equation with the 
/^-function for the invariant charge. 

3 This choice is the following: 7 = 70 = — 1; 7$ = 0, i > 1, where 7« is the coefficient by the 
(i + l)th power of the perturbative running coupling on the right-hand side of Eq. ([l]). 
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V 1 ) 

an 



(Q 2 



An z - 1 

(3q z In z' 



A 2 ' 



(7) 



where /?o = 11 — 2n//3 is the first coefficient of the /3-function. At the higher loop levels 
there is only the integral representation for N a an (Q 2 ). So, at the i-loop level we have 



an 



(Qo)— exp 



(7 



111 



(T 



da 



(8) 



ai*'(- a — is) 



/(2m) is the spectral density, and z 



where p®(cr) '- 
Ql/A 2 is the normalization point. 

Figure [1] shows the new analytic running coupling computed at the one-, two-, and three- 
loop levels. It is clear from this figure that our analytic running coupling possesses the higher 
loop stability. Moreover, it can be shown that the singularity of the new analytic running 
coupling at the point Q 2 = is of the universal type at any loop level. This is clear from 
the following simple consideration. When z — > the basic contribution into Eq. (|j) affords 
the integration over the small a region. The spectral density p^ i '(u) at any loop level has 
the same limit when a — > 0: p^(er) — > p^(cr) = (In 2 a + n 2 )^ 1 ||. Hence the new analytic 
running coupling (§) has the unique behavior when z — > 0. 



III. QUARK- ANTIQUARK POTENTIAL GENERATED BY THE NEW 
ANALYTIC RUNNING COUPLING 

Here we are going to use the new analytic running coupling for obtaining the interquark 
potential. We proceed from the standard expression ]5|,|TT]] for the qq potential in terms of 
the running coupling a(q 2 ), 

. 167T f°° a(q 2 ) e iqr , 
V(r) = -— / -^tjtt^ dq. 9 
3 Jo q l (2ix) s 

For the construction of the new interquark potential V(r) we shall use the new analytic 
running coupling (0) 

p zmz A 2 

Upon the integration over the angular variables and the substitution q/A — > q, rA — > R in 
Eq. @, one gets 

32 ~ 

N V(r)=-—A.V(R), R = Ar } (11) 

where 

t~,/t^ f 00 o 2 — 1 sin(gi?) , , . 

V ( R ) = 21 2 p 1 dq 12 
Jo In q 2 qR 

is the dimensionless potential. 
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In order to perform the integration in Eq. (|T2| ) we consider the auxiliary function 



I(n, R) = lim / 



Q 



a^o+ Jo ln(a + q 2 



sm(qR) dq. 



(13) 



Here the parameter a is introduced for shifting the origin of the cut along the imaginary 
axis Im q. It is obvious that 



V{R) 

For even n the integrand in Eq. (| 



1 

7? 



[I(0,R)-I(-2,R)}. 



is an even function of q. Therefore 



(14) 



I(n, R) = — Im lirn^ J(n, a, R), 



where 



/oo 
F{q) dq, F{q) 
-oo 



n n—l„iqR 

ln(a + q 2 ) 



(15) 



(16) 



The sign V means the principal value of the integral. 

The function F(q) in Eq. (ITBl) has the cuts (—zoo, —it/a], [iy/a, zoo) and simple poles at 
the points q = =F\A ~~ a - Let us consider the integral of the function F along the contour 
T shown in Fig. ||. The function F(q) has no singularities inside the contour T, therefore 
fpF dq = 0. Contribution to this integral of the semicircle of infinitely large radius in upper 
half-plane (see Fig. 0) vanishes. Performing the integration along the two semicircles c_ and 
c + of the vanishing radius and along the cut C on the imaginary axis, we obtain 



J{n, a, R) 



%TX< 



+2i 



71-2 



dx 



\*fa ln 2 (x 2 — a) + 7T 2 j 

Hence, for even n the function J(n, R) in Eq. flT3]) takes the form 



I{n, R) 



71 



— cos R 



■l) n/2 M{R,n) 



where 



M{R,n) 



o In" X 2 ) + 7T 2 



dx. 



(17) 



(18) 



(19) 



It is rather complicated to perform the integration in Eq. (|19|) explicitly. Therefore we 
address the study of the asymptotics. First of all, we would like to know whether the qq 
potential N V(r) in Eq. (O) provides the quark confinement. For the investigation of the 
potential behavior at large distances it is enough to consider the asymptotic of the function 
J\f(R, n) in Eq. (|T9|) when R — > oo. This function can be represented in the following way 
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gn , 



-Rx 



dR n Jo x 



\I1 2 (X 2 ) + 7T 2 



dx. 



(20) 



At large R the basic contribution into Eq. ( ^Op gives the integration over the small x region. 
Let us transform M(R, n) identically: 



Af(R, n) 



dR n Jo 4x(\n 2 x + it 



l) n e~ m \ l 3L 



1 + L 



dx. 



(21) 



where L = 7r 2 /(41n 2 x). Neglecting the second term in the square brackets in Eq. fl21"D , we 
use the formula (4.361.2) from Ref. 



o x(ln X + 7T 2 



-dx = e M — ^(/x), Re /i > 0, 



(22) 



where z/(/z) is the so-called transcendental //-function [jl6 



Z/(/i) 

Eventually, we obtain for R — > oo, 

(-l) n 



o r(t + i)' 



A/"(.R, n) 



e R - u(R) + 



o r(t + i) 



(23) 



(24) 



Taking into account Eqs. (|14l), (0), and (f?|) one can present the quark-antiquark po- 
tential (|TT|) at large i? in the following way: 



A? 



V{r) 



8tt A 
3Po~R 



2 (e fl - u(R)) + J 



2 i2*dt 

nt+i) 



(25) 



The behavior of the potential V(r) at r — > oo is determined by the last term in Eq. (p5|).f| 
Integration of this term by parts gives 



2 R l dt 



R 2 



o r(t+l) 21ni? 

1 

_ m~R 



1 + 2VM 



i + E 



/i(0) 



(26) 



where 



fn(t) 



d n (-l) n 

ds n r(s + 1) 



(27) 



4 It follows directly from the asymptotic of v{R) (see Ref. J16||), and from a simple reasoning. 
Really, if R > the term f(R) = e R — v(R) is non-negative and f'(R) < 0. Hence, /(-R) — ► const 
when R — > oo, and its contribution to ^V(r) at large i? is of l/i?-order. 
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In the limit R — > oo, Eq. (p6|) takes the form 

R?dt R 2 



T(t + 1) 2\nR' 

Therefore the quark-antiquark potential N V(r) proves to be rising at large distances 



(28) 



iVy( r )^^l A .I_^_ r^oo, (29) 
V ; 3/5 21n(Ar)' v ' 

Thus the new analytic running coupling N a an (g 2 ) [see Eq. (|7|)] leads to the rising quark- 
antiquark potential N V(r) which can, in principle, describe the quark confinement. 

It is important to point out that the behavior of the potential N V(r) when r — > has 
the standard form determined by the asymptotic freedom (see, e.g., Ref. [11]]), 



^W-^A-— i— , r^O. (30) 

Unfortunately, it is impossible to obtain the explicit dependence V(r) for the whole region 
< r < oo. A simple interpolating formula, which can be applied for practical use, will be 
given in the next section. 



IV. DISCUSSION 

Let us discuss briefly the higher loop contribution. As was mentioned in the Sec. II, the 
singularity of i-loop analytic running coupling N a^(q 2 ) at the point q = is of the universal 
type at any loop level. Therefore, when q — > we have N a^(q 2 ) ~ N a^(q 2 )C\ where 
C l are constants. Taking into account that the maximal difference between N a^(q 2 ) and 
N a^(q 2 ) is in the small q 2 region, we arrive at the following conclusion. The account of the 
higher loop corrections leads to changing the slope of the qq potential N V(r) when r — > oo. 
This corresponds to a simple redefinition of the parameter A in Eq. (|29|) at the higher loop 
levels. 

As far as the scheme dependence of this approach, we have to point out the following. 
It was shown in || that the solutions of the analytized RG equation at the higher loop level 
have extremely weak scheme dependence. In particular, the solutions of the RG equation 
with MS and MS schemes, are practically coinciding. Hence, at the higher loop level (there 
is no scheme dependence at the one-loop level), the use of different subtraction schemes 
leads to the slight variation of the qq potential. 

Thus, neither higher loop corrections, nor scheme dependence, can affect qualitatively 
the result obtained in the previous section. 

For the practical use of the new potential it is worth obtaining a simple explicit expres- 
sion that approximates it sufficiently well. For this purpose one can use, for instance, the 
approximating function 




\n 2 R 




(31) 
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which has no any unphysical singularities and possesses the asymptotics ( p9|) and (p0|). This 
function is obtained by smooth sewing the asymptotics 



S7T 



N V(r) ~ --A 



57T 



^y(r) ~ —A 

1 ; 3/3 flln(i2) 



2\n{R) \n 2 (R)_ 
1 



R 



oo, 



i? _> 0, R = Ar. 



(32) 
(33) 



The formula (]32f) keeps explicitly the second leading term of the expansion fl26|) , /i(2) = 
0.461. Some terms have been introduced into Eq. (|31|) only for eliminating the singularity 
at the point R — 1. It should be mentioned here that the next terms in the expansion ( |26|) 
practically do not affect the shape of U(r). Of course, the function ( |3T| ) is not the unique 
interpolating function between asymptotics (j32j) and (|33|). Nevertheless, the comparison of 
U (r) with the phenomenological potential 



c 



... , 4a 

1/ (r) = h err + const 

3 r 



(34) 



(the so-called Cornell potential M) shows their almost complete coincidence (see Fig. [|). 
The fit has been performed with the use of the least square method in the physical meaning 
region 0.1 < r < 1.0 fm [||]. The varied parameter in Eq. (|3l|) is A. The possibility of 
shifting the potential V(r) in Eq. ( ^4[) by a constant was also used. A rough estimation 
of A in the course of this fitting gives A ~ 500 MeV. This is in agreement with the values 
obtained earlier in the framework of the analytic approach to QCD [H]. 



V. CONCLUSION 

In the paper the quark-antiquark potential is constructed by making use of the new 
analytic running coupling in QCD. This running coupling arises under analytization of 
the renormalization group equation before its solving. The rising behavior of the quark- 
antiquark potential at large distances, which provides the quark confinement, is shown 
explicitly. The key property of the new analytic running coupling, leading to the confining 
potential, is its infrared singularity at the point q 2 = 0. At small distances, the standard 
behavior of the potential, originating in the QCD asymptotic freedom, is revealed. It is 
also demonstrated that neither higher loop corrections, nor scheme dependence, can affect 
qualitatively the obtained result. The estimation of the parameter A in this approach gives 
a reasonable value, A ~ 500 MeV. 

In further studies it would undoubtedly be interesting to consider in this approach the 
dependence of the qq potential on the quark masses. 
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FIG. 1. The normalized new analytic running coupling R(z) = N 5 an ((5 2 )/ N 5 an ((5o) at the 
one-, two-, and three-loop levels. The normalization point is Q% = 10A 2 , z = Q 2 /A 2 . 
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FIG. 3. Comparison of the potential U(r) given by Eq. ( |3l|) (solid curve) with the phenomeno- 
logical Cornell potential (o), Eq. (|34j). The values of the parameters are: a = 0.39, a = 0.182 GeV 2 
§], A = 530MeV, n/ = 5. 
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